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Abstract 

We estimate the ratio R = (73/(72 of the critical coupling constants g 2 and 53 which 
are required to achieve binding of 2 or 3 bosons, respectively, with a short-range 
interaction, and examine how this ratio depends on the shape of the potential. Simple 
monotonous potentials give R ~ 0.8. A wide repulsive core pushes this ratio close 
to R = 1. On the other hand, for an attractive well protected by an external repulsive 
barrier, the ratio approaches the rigorous lower bound R = 2/3. We also present 
results for N = 4 bosons, sketch the extension to N > 4, and discuss various 
consequences. 



I. INTRODUCTION 

The phenomenon of "Borromean" binding is well-known [jl]J^]. In our world with 3 dimen- 
sions, a short-range potential gv (r) does not always achieve binding of two bodies, even if v is 
attractive or contains attractive parts. A minimal strength g 2 is needed, more precisely mg > g 2 , 
where m is the mass of the constituents. Similarly, binding three identical bosons of mass m re- 
quires mg > #3 for the pairwise interaction g v ( r 'ij), where denotes the distance between 
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particles % and j. The crucial observation is that g 3 < g 2 , implying that for a coupling g such that 
93 < 9 < 92, the 3-body system is bound while none of its subsystems is bound. 

An example is the 6 He nucleus, considered schematically as a (a, n, n) system. It is sta- 
ble against spontaneous dissociation, while neither the (a, n) = 5 He nor the (n, n) systems are 
bound. The name "Borromean" was given to such nuclei after the Borromean rings, which are 
interlaced is such a subtle topological way, that when one removes one of them, the two others 
become unlocked. 

Borromean binding is implicit to understand the Thomas collapse [[5fl. When the range of 
the potential v is reduced, or equivalently, when g — > g 2 from above, the 3-body binding energy 
E 3 (g) becomes very large compared to the 2-body energy E 2 . Also the Efimov effect [Q], i.e., the 
proliferation of loosely-bound excited states in the 3-body spectrum near g = g 2 implies that the 
3-body ground-state already exists at this point. 

Let R = (73/(72 be the ratio of critical coupling constants. For simple monotonous potentials, 
such as Yukawa, Gaussian or exponential, it is found [|5]] that R is very close to 0.8. This is in 
agreement with the rigorous lower bound R > 2/3 [Q]. The fact that all simple potentials give 
almost the same R ~ 0.8 is understood as follows: at vanishing energy, the wave function extends 
very far outside the potential well, and thus does not probe very accurately the details of the short- 
range interaction, which is just seen as a contact attraction. 

There are, however, reasons to believe that R can appreciably differ from 0.8. The aim of 
the present paper is precisely to study how R evolves when one starts from a simple monotonous 
potential and adds either an inner core or an external barrier. 

When an external barrier of growing size is added to the potential, R evolves from R ~ 0.8 to 
R — > 2/3. An example is provided by v oc r 2 exp(— 2/ir) — exp(— /xr) when /1 varies, or, similarly, 
by combinations of Gaussians. 

When an inner core is implemented, a transition is observed from R ~ 0.8 to R — > 1. This will 
be seen for the Morse and the Poschl-Teller potentials when an appropriate parameter is varied. 
An extreme case consists of an hard core of radius c and an attractive delta-shell — S(r — d) located 
at d > c. The critical strength g 2 can be calculated exactly. One can also calculate exactly the 
strength g^ which makes the 2-body scattering length vanish and hence is sufficient to bind the 
infinite boson matter [Q], with the result g^j g 2 = c/d. Thus, as d — > c, g^j g 2 approaches 1 and 
so does any g N / ' g 2 ratio with iV finite. 

Note that the ratio R cannot exceed R = 1 for the additive potential V=g ^ vfaj), provided 
v is purely attractive. This means that one cannot conceive a situation where 2-body systems are 
bound while a 3-boson system is unbound. The following proof is due to Basdevant [|7p. For 
g < g 2 , let (p(r) be the ground state wave-function of the 2-body system, with energy E 2 . The trial 
wave-function \I> = V 9 ( r i2)v 9 ( r i3) can be used for the 3-body Hamiltonian written as, 

H=4(- + -)+4(- + -) + — + + + v(r 31 ), (1.1) 
2 \mi m 2 J I \vtlx m 3 J mi 

leading to an expectation value 2E 2 if the interaction v(r 23 ) is neglected. So E 3 < 2E 2 < if 
v < 0, q.e.d. The proof holds for an asymmetric interaction V=J^ i>y- with v\ 2 binding (1, 2) and 
fi 3 binding (1, 3), and v 23 only weakly attractive or even vanishing. We believe that this result 
remains true if v is not purely attractive, but we have not been able to prove this generalization. 

This paper is organized as follows. In Sec. we discuss how to compute accurately the 
critical couplings g^. In Sec. [n|, we present the results obtained for the Morse model and a 
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few other potentials: the three-body binding energy obtained at g = g 2 , at the edge of binding 
for the two-body systems; the critical coupling g 3 for three-body systems; an estimate of g A , the 
minimal strength necessary to bind four bosons. It is expected that g 4 < g 3 , with, however, the 
constraints (74/(73 > 3/4 and (74/(73 > 1/2 established in Ref. [Q]. The numerical estimate of #4 
requires delicate variational calculations, especially when the potential displays both attractive and 
repulsive parts. A simple extrapolation to g^, i.e., the infinite boson matter, will be presented in 
Sec. [IV]. Some conclusions and a list of open problems are presented in Sec. |VJ 

II. VARIATIONAL METHODS 

There are well-known techniques, in particular variational methods [||], to compute with very 
high precision the binding energy E N (g) of a system of iV particles in a regime g > g N where 
binding is established. It is a slightly different art, however, to estimate the value g N corresponding 
to the border of the stability domain. Even in the simple case of N = 2 constituents, this is not 
completely obvious, as seen, e.g., from the discussion in Refs. [p[,|IO|] for the Yukawa potential. 

A first strategy consists in computing accurately the binding E N (g) in a domain where binding 
occurs and letting g decrease. As a behavior 

En oc —(g — gN) 2 (2.1) 

is expected, one better looks at (— E N ) 1 ^ 2 as a function of g and checks a straight behavior as 
E N — > 0. As in Ref. [JTlQ, a Pade-type of approximation is found adequate to extrapolate towards 
gN- In a typical variational method, the Schrodinger equation (T + V)^ = En 1 ^ is solved by 
expanding the wave function on a basis of functions 

i 

In a given set of cfi, the weights C-i (represented by a vector C) and the variational energy E are 
obtained from a generalized eigenvalue equation 

[f + gV)C = ENC, (2.3) 

involving the restrictions of the kinetic energy operator T and potential energy V to the space 
spanned by the ipi, and a definite-positive matrix N, which does not reduce to the unit matrix 
when the cpi are not orthogonal. 

An alternative (though not strictly legal) method for estimating the critical coupling gN consists 
of looking directly at the point E = and rewrite the eigenvalue equation as 

VC = — TC, (2.4) 

9 

involving, again, an Hermitian matrix on the left, and a definite-positive matrix on the right. In 
principle the wave-function needs not be normalizable at E = 0, but in practice, one can use a 
basis of normalizable functions, provided one allows for components with very long range. 

The results presented below have been checked using both the extrapolation method and the 
direct estimate of the critical coupling. 
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When the number of terms in Eq. Q2.2| ) is incremented, there is a dramatic increase of the 
number of non-linear parameters entering the basis functions (the coefficients in the examples 
below). The minimization of the variational energy by varying these parameters becomes i) am- 
biguous, as neighboring sets of values give comparable energies, and ii) intractable, even with 
sophisticated minimization routines. A simple trick [|T2p is inspired by the work of Kamimura 
[|T3[]. It consists of imposing all parameters to be chosen in a single geometric series. Then 
only the lowest and the largest values have to be optimized numerically. The minimization is much 
faster. The slight loss in accuracy is more than compensated by the possibility of increasing easily 
the number of terms. This works rather well for achieving a reasonable accuracy. When one aims 
at very precise results, more sophisticated techniques are required, such as the well-documented 
and powerful stochastic variational method (SVM) [^i[l4|]. 

It remains to choose the basis functions in Eq. (|2.2|). We have compared the results obtained 
with exponential functions 

(pi = exp f - ^ a ij r ij) H ( 2 - 5 ) 

and Gaussians 

ifi = exp (- a ij r l) H > ( 2 - 6 ) 

where the parenthesis can be rewritten as the most general quadratic form involving relative Jacobi 
coordinates. In both cases, the dots are meant for terms deduced by permutation, to ensure the 
proper symmetry properties of the trial basis. 

The former basis is by far more efficient when the expansion fl2.2|) is limited to a small number 



of terms. For instance, a single exponential function is sufficient to demonstrate the stability of 
the ion Ps = (e + e~e~) in quantum chemistry, while several Gaussians are needed. However, 
when the number of terms increases, the exponential basis, even when associated with a stochastic 
search, tends to give rise to numerical instabilities similar to those described, e.g., by Spruch and 
Delves QT3|] . The problem can certainly be circumvented [|T6|], but we found it more convenient to 



use SVM with Gaussians to get stable and accurate results. Anyhow, the results involving more 
than N = 3 particles have been obtained with Gaussians only, since one cannot derive simple 
analytic expressions for the matrix elements within the exponential basis. Note that when N 
increases, the surface and tail of the system play a relatively less important role, so the use of 
Gaussian functions should become more appropriate. 

IE. RESULTS FOR g 3 /g 2 AND g 4 /g 2 

In this section, we present some results on R = g^/gi and i? 4 = (74/(72- We restrict ourselves to 
symmetric 3- or 4-body systems, involving identical bosons. Some results on equal-mass particles 
with asymmetric interaction have been given in Refs. [p|,^7[]. 



A. Monotonous potentials 

We consider here three simple functional forms, Yukawa, exponential and Gaussian, corre- 
sponding to 
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— v(r) = , expf — fj,r), exp(— /i r ), (3.1) 

r 

respectively, where, without loss of generality, the range parameter fi and the constituent mass m 
can be set to /z = m = 1 by simple rescaling. 

The critical couplings are displayed in Table |. As already stressed, the most remarkable feature 
is the close clustering of all values of R near 0.80. This means a 20% window for Borromean 3- 
body binding. Similarly, all values of R± are found around 0.64. There is no obvious meaning to 
the observation that (7.4/(73 ~ (73/(72- Anyhow, gN+i/gN cannot be smaller than N/ (N + 1) [Q], so 
it should tend to 1 as N increases. 



B. Potentials with external barrier 

The potential 

v(r) = aexp(-//V/2) + 6exp(-2/iV) (3.2) 



has been used by Nielsen et al. pl[|, to study Borromean binding in two dimensions. By rescaling, 
one can fix m = \i = 1. The cases (a, b) = (—1, 0) and (a, b) = (1, —2) are shown in Fig. [j] 
for illustration. For a and b both negative, this potential reduces to a simple monotonous function, 
and, not surprisingly, a ratio R ~ 0.80 is obtained for the critical couplings g% and gi. With a and 
b of different signs, one can build a potential which looks like an almost pure harmonic oscillator 
at small values of r and vanishes only at distances which are very large as compared to the size of 
the ground state wave function. One then obtains R — > 2/3. 

In this limit, we are, indeed, approaching the situation where the decomposition 

H 3 (m,g) = J2H?' j \3m/2,g) (3.3) 

corresponds to an exact factorization of the wave function [|l9i^tl, and thus the vanishing of 
H 3 (m,g) implies that of H^^ {2>m/2, g) i.e., g s = 2g 2 /3 [j5J]. Otherwise, one simply gets 
from Eq. (|3~3[), when saturated with the exact 3-body wave function, E 3 (m,g) = ^2(H 2 ) > 
3£ 2 (3m/2, (?), i.e., g 3 > 2g 2 /3. Here 

(3.4) 

k k<l 

is the translation-invariant Hamiltonian describing the relative motion of N particles. 

Table | gives the values for a pure Gaussian potential, corresponding to (a, b) = (—1, 0): one 
obtains g 2 = 2.680, (73/(72 = 0.79 and (74/(73 = 0.80. For the potential (a, b) = (1, —2) also shown 
in Fig. [1|, the values become: g 2 = 21.20, (73/(72 = 0.672 and (74/(73 = 0.754. We are already very 
close to the limit where gN+i/gN — N/(N + 1). 



C. Morse potential 

The Morse potential reads 
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v(r) = exp [— 2/i(r — r )] — 2 exp [— /x(r — r )] 



(3.5) 



Again, one can set m = /x = 1 by rescaling. The shape is displayed in Fig. |2|, for r = 1 and 
r = 2. The 2-body problem with this potential can be worked out exactly [[TSp. In particular, the 
critical coupling g 2 is obtained from an equation involving the Kummer function, which can be 
solved easily. Our normalization is such that g 2 — > 1/4 as r increases. The critical couplings g 3 
and <?4 have been estimated numerically, as well as E 3 (g 2 ), the 3-body ground state energy at the 
edge of binding 2-body systems, and similarly E^(g 3 ). The results are shown in Table pi 

A warning is that the calculation becomes very difficult as r becomes larger than about 3. Our 
parametrization becomes inadequate. The vanishing of the wave function at small interparticle 
distance r\j is obtained at the expense of huge cancellations in the expansion fl2.2p. This consider- 



ably reduces the accuracy. Specific methods can be developed for interactions with hard core, see, 
e.g., and references therein. Our results, however, seem good enough to show unambiguously 
the trends of the Rn ratios as the size r of the core increases. 

The size of the attractive pocket is measured by the interval between r\, where the potential 
vanishes, and r , where it reaches its minimum. Within our normalization, Sr = r — r\ = In 2 
is constant. As r increases, Sr/r — > 0, and the Morse potential becomes similar to the attractive 
delta-shell with hard core described in Sec. | in the limit d/c — > 1, and then a behavior R — > 1 is 
expected. This is clearly observed in Table 0. The trend is, however, rather slow. For moderate 
values of r , our numerical results are well reproduced by a fit 

Sr 0.43 

R ~ 1 - c— ~ 1 , (3.6) 

with c ~ 0.62. For larger r , some departure is observed, probably due to the difficulties in the 



variational calculation. We believe the behavior (3.6) is rather general 



D. Poschl-Teller potential 

It reads 

= q(q-l) _ q(q + l) ^ ^ ?) 

sinh 2 (/ir ) cosh 2 (fir) ' 

with, again, \i — m—\ for the range and the mass of each constituent. The strength factors of the 
repulsive and attractive terms are tuned to give a zero-energy 2-body state, i.e., g 2 — \ [[[Sp. The 
potential is drawn in Fig. [3] for some values of a. 

In the case a = 1, we have a simple monotonous potential, and, not surprisingly, a value 
R ~ 0.8 is found, as seen in Table p|. 

For a > 1, we can again define the size of the attractive well as Sr = r — r 1? with v(ri) = 



and v'(r ) = 0. The fit of R using the empirical formula ( pT6[ ) turns out to be quite good. This can 
be checked from the values displayed in Table |I| 

For both the Morse and Poschl-Teller potentials, the approximate equality of (74/^3 and g^/ g 2 
survives a strong repulsive core, unlike the case of an external barrier. There is, however, a slight 
difference in the patterns exhibited by these potentials. In the Morse case, the ratios g^/g^ and 
93/92 start departing from about 0.80 at the same value of the parameter r for which g^j g 2 
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becomes positive, i.e., binding the infinite boson matter requires a minimal strength. In the Poschl- 
Teller case, (74/(73 and (73/(72 immediately increase when the parameter a becomes larger than 1, 
though (7oo / 'g 2 still vanishes for a while. 



E. More complicated potentials 

In Ref. [O] the ratio R is studied for the potential 

v(r) = r 2 exp(— 2/xr) — exp(— fir), (3.8) 

which is shown in Fig. |], for selected values of the range parameter fx. For fx very large, it is almost 
attractive, and one thus obtains the usual R ~ 0.80. On the other hand, for very small fx, we have 
an almost pure oscillator v(r) ~ r 2 — 1 in the region of interest, and one gets R — > 2/3. For 
intermediate value of fx, the potential exhibits both an internal pocket of attraction and an external 
attractive tail. 

The numerical results are shown in Fig. |5|. The expected behavior i? = 2/3for/i^0 and 
R ~ 0.80 for fx — > 00 are verified. Near fx = 0.70, there is an interesting tunnelling between the 
internal and the external pockets of attraction. The barrier is seen as an internal core by the latter, 
and this pushes R toward 1, as for a Morse potential of large radius. 



IV. LARGER SYSTEMS 

The case of an infinite boson matter sheds some light on our discussion. For a purely attractive 
potential, a system containing many bosons is bound, however weak is the strength g: we will thus 
set (7oo = 0. Now, if the potential contains a large repulsive part, it is conceivable that binding 
requires a minimal strength of the potential, say g > > to pull the wave function in the 
attractive parts of the potentials. A result which looks a posteriori reasonable, is that g^ is the 
value of the coupling for which the scattering length vanishes [j^]. Indeed, the optimal state of the 
infinite boson matter is a compromise between the large-density limit, for which the kinetic energy 
is too large, and the extreme dilution. The latter case, dominated by two-body collisions at zero 
energy, should exhibit a tendency towards binding, i.e., a negative scattering length. 

The scattering length can be calculated analytically or numerically for the potentials considered 
previously. Then it is rather straightforward to determine the value g^ of the strength which makes 
it vanish. Of course, for a potential whose integral is negative, the scattering length is already 
negative in the weak coupling limit, and remains negative until g = g 2 . This corresponds to 
9oo = 0. 

In Tables and[ni|, we display the value g^ for which the scattering length vanishes. This is 
simply #00 = for the monotonous potentials of Table | and the limiting cases r = of the Morse 
potential and a = 1 of Poschl-Teller. As r or a increases, one observes almost simultaneously 
g^ becoming finite and R = #3/(72 departing from about 0.80 and approaching R — 1. This is the 
limit between, say, simple potentials which are purely attractive or contain a small repulsion, and 
non-trivial potentials with a strong core. 
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V. DISCUSSION AND OUTLOOK 



In this paper, we have studied some aspects of the phenomenon of Borromean binding in three 
dimensions by comparing the critical couplings g N required for binding iV = 2, 3 or more bosons 
interacting with various types of potentials. 

All monotonous, short-range potentials give almost the same ratios gj,j g 2 — 0.80 and g±/gz — 
0.64. 

We then considered potentials with a short-range attraction and an external repulsive barrier, 
which behave very much like an oscillator, and, not surprisingly, give ratios of critical couplings 
close to the lower bound 9n+i/9n = N/ (N + 1). These potentials with an external barrier are 
not very often encountered in physical systems. They are, however, interesting, since, according 
to Ref. pi"|], they are the only ones to give rise to Borromean binding in two dimensions. 

We then studied the more physical case of potentials with a strong repulsive core at short 
distances. The window for Borromean binding turns out to be much narrower than for purely 
attractive potentials. 

The present investigation could be extended to excited states. In particular, as long as g 3 < g 2 , 
the Efimov effect should remain as g approaches g 2 . It would be interesting to study how the onset 
and disappearance of Efimov states change when the strength of the core is varied. 
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TABLES 



TABLE I. Comparison of critical couplings gpj to achieve binding of N = 2, 3 and 4 identical bosons 
in a Yukawa (Y), exponential (E) or Gaussian (G) potential. Also shown is which correspond to a 
vanishing of the 2-body scattering length, and the 3 -body energy E% obtained for g = g2, and the 4-body 
energy E 4 for g = g 3 . 
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TABLE II. Same as Table |[ but for a Morse potential (^5), whose characteristic radius tq is varied 
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TABLE III. Same as Table |, but for a Poschl-Teller potential ( p.7[ ), for several values of its parameter a. 
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FIG. 3. Shape of the Poschl-Teller potential for a = 1 (monotonous) and, from left to right a = 2, 3, 5 
and 9. 




FIG. 4. Potential of Eq. (^8|), for n = 0.4 (thick line), /u = 1 (thin line) and fj, = 0.65 (dotted line). 
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FIG. 5. Computed value of R = g^/gi for the potential of Eq. ( |3.8[ ), for various values of the range 
parameter /i. 
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